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Abstract. We propose a gauge-independent mechanism for the area-law behavior of Wilson loop expecta- 
tion values in terms of worldsheets spanning Wilson loops interacting with the spin foams that contribute to 
J^j^' the vacuum partition function. The method uses an exact transformation of lattice-regularized Yang-Mills 

i—{ ' theory that is valid for all couplings. Within this framework, some natural conjectures can be made as to 

what physical mechanism enforces the confinement property in the continuum (weak coupling) limit. Details 
for the 5(7(2) case in three dimensions are provided in a companion paper. 
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1. Background 



Our current understanding of confinenient in Yang-Mills gauge theories encompasses a seemingly diverse 
, collection of approaches, including center vortices [37|, dual superconductors [28 l f33 l f38] . and effective strings 

of various types [Sl[32] — to name just a few that have attracted interest over the years. 

In the present work we propose a mechanism for confinement that is similar in concept to the well-known 
(-H ' proof of confinement in the (unphysical) strong coupling limit [30j . but suitably generalized so that the 

analysis is valid at the arbitrarily weak couplings that are characteristic of the continuum limit. Our main 
result is a new form of the Wilson loop criteria for confinement in lattice gauge theory, in terms of which we 
^ , conjecture a general confinement scenario that can be applied to any specific Yang-Mills theory (i.e. choice 

Q>^ ' of gauge group and space-time dimensionality). 

00 . We refer the reader to a companion paper [llj for a demonstration of the framework developed herein 

OO ■ for the case of G = SU{2) in three dimensions. The present paper focuses on establishing the worldshcet- 

background decomposition of spin foams and outlining in general terms those aspects of our proposal that 
OO ' are model independent. For concreteness, we do however specialize to the case of SU{2) starting in Section 

, 3. The scenario is also well-defined in the SU{3) case, but more detailed analysis will require further 

developments in the evaluation and asymptotic limits of SU{3) spin networks. 

The structure of the present paper is as follows. In Section 2, we shall review the standard lattice 
regularization of pure Yang-Mills theory on the lattice, the definition of a Wilson loop observable, and the 
concept of character expansion for the conventional amplitude. This will allow us to give a brief account of 
^ ' the widely known proof of confinement at strong coupling and fix basic notation to be used throughout. 

5^ , In Section 3, we will introduce the spin foam formulation of lattice Yang-Mills theory, which is exactly 

dual to the conventional formulation. By analysis of admissibility conditions, we show that contributions 
to the Wilson loop expectation value can be organized as a double sum over worldsheets connected to the 
Wilson loop and over vacuum spin foams, a result already implicit but only realized in a limited way in past 
work [23; on strong coupling. In Section 4 we show how this decomposition of spin foam configurations can 
lead to a natural generalization of the strong coupling analysis of confinement. We also discuss how certain 
features of a roughening transition are present in transitioning to the weak coupling limit but present no 
obstacle within the spin foam framework, unlike past methods that use strong coupling expansions. We 
conclude in Section 5 with a discussion of our results and the further lines of inquiry that they suggest. 

2. Review of lattice formulation and confinement at strong coupling 

Recall the Euclidean partition function of pure (non-abelian) Yang-Mills theory in D dimensions, with 
gauge group G = SU{N) for N > 2. In the continuum it takes the form 

(1) Zv^ jvA oM-S[A]), 



where A = A'^{x) is the gauge field, S is the action functional, and VA is the fimctional integration measure. 

Definition 1 (Lattice). We shall denote our oriented D— dimensional hypercubic lattice by k, and its 
associated set of vertices, edges, and plaquettes by V, E, and P, respectively. 

We apply the standard lattice regularization by defining a functional integral measure as an integral over 
a product of \E\ copies of G using the Haar measure: 

(2) - n ^5e- 

eeE 

To proceed we replace the continuum action functional by a discretized action, S = S[g] that reproduces the 
continuum action S[A] in the limit of zero lattice spacing. A discretized action functional S can be formed 
by summing over a local function S{gp) of variables at a plaquette, S[g] = J2p£P ^idp)^ where the group 
element gp is the holonomy around an oriented fundamental plaquette p. That is, gp = 5i52ff354j where gi is 
either the group element assigned to the ith edge of p or its inverse if the orientations of p and the ith edge 
are opposing. This yields the conventional lattice partition function 



(3) Zv= (ll dg. 



\eeE 



Let F denote a lattice loop (a closed, non-self-intersecting lattice curve). We define a Wilson loop to be a 
pair (F, w) where F is a lattice loop and w is the label of the irreducible representation in which the trace 
is to be taken. For brevity we denote this data by F^, = (F,-u;); the object F„ can be thought of as a loop 
in the lattice colored by the charge w. The expectation value (Or„) associated with F^, is defined as a field 
insertion into the partition function as follows: 

,^ , _ !{Y{eeEd9e) Tr(n,gr^^(ge))e-^p^^^(^^) 

Here Jleer ^wige) denotes a product of representation matrices D^^g) in the wth irreducible representation 
of G (in studying confinement for G = SU (2), the case w = ^ is usually considered). We can now give a brief 
sketch of the proof of confinement (at strong coupling) by way of the character expansion. This approach 
to expanding the partition function is an exact transformation, as the summation over all diagrams yields 
a convergent sum for any non-zero value of the coupling constant j30| . We have for the vacuum partition 
function the following expansion: 



(5) Zv = 



/fn^5e)n(E^>x>(5p)| 

\e&E ) pep J 

fl f [U'^Se] l[c,^xjgp), 

b„) Vee-E / PGP 



{jp} \ee-E / pel 

where e~^^Sp'> — '^jpXjpidp) defines the character coefficients cj^ associated with the chosen plaquette 

action S{gp). In the second line we have interchanged the order of summation and integration operations. The 
notation {jp} refers to the set of all labellings of the plaquettes of k by unitary irreducible representations 
(irreps) of G. Each such labelling can be viewed as a "diagram" — however, only a special subset of 
diagrams in the summation are non-zero. The labelling^ which give non-zero (we refer to these as admissible) 
contributions to Zy satisfy certain constraints on the representations of plaquettes that meet at an edge. 
In SU{2), a geometric characterization is possible: admissible diagrams can be viewed as collections of 
surfaces that are closed (no free edges allowed), branched (can meet other surfaces along vertices, edges, and 
plaquettes) and colored by the group irrep labels [TBlfTTlIM] . 



Note we are not following the more traditional method of expanding in powers of /3, but rather products of character 
coefficients (see for example I17| for more on this approach to strong coupling expansions). 
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By applying the same procedure of character expansion described for the vacuum partition function above, 
one can associate a charged partition function with the numerator of the expectation value (|4]) : 

(6) Zr ^ |(^n^-9^^) Tr(^jy;i?.(.9,)^e-^.ep^(^'') 

= E/(n'^5e) Tr(ni^i(.9e)) n^^^^.M- 

{jp} VeGE / Veer /pGP 

As a function of Wilson loop F, the numerator of the expectation value (Or„) can thus be re-expressed as 
a sum over diagrams labelled by the unitary irreducible representations of G as: 

(7) (OrJ - 1^ - ^y' E / ( n <9e ) Tr ( [] ^^(.9-) ) H ^..^.(Sp). 

The summation is again over all labellings of the plaquettes of the lattice by unitary irreps of G. However, due 
to the presence of representation matrices associated with the Wilson loop, the admissibility constraints are 
locally modified along F. One important consequence of this is that for any contribution to the numerator of 
(Or„ ) there is a lower bound on the number of plaquettes carrying non-trivial irrep labels that is proportional 
to the area of a minimal surface spanning the loop. As we shall argue more rigorously below, the expectation 
value of a Wilson loop operator can be expressed as a ratio of two partition functions, involving two different 
collections of surfaces: one in which all surfaces close (Zy) and a second in which a surface always ends on 
the Wilson loop, while all other surfaces close (Zy). 

Transition to the continuum limit is controlled by /?, a parameter inversely proportional to the coupling 
constant. Although suppressed in the above notation, /3 enters through the character coefhcients Cj^ = Cj^ (/?). 
Later, we shall see how the character coefficients become plaquette amplitudes in the spin foam formulation, 
and the fact that they contain all the /9-dependence of the model will be important to understanding the 
weak coupling limit. Returning now to ([7|), it can be shown that for sufficiently small /3, the strong decay of 
for large spin combined with bounds on the integrals of products of characters X] allows for a series 
expansion for (Or„) having a finite radius of convergence and leading order terms that are exponentially 
decaying in the area of the Wilson loop [S^ . 

As is typically remarked when considering the extension of these arguments to weak coupling, increasing 
/? (as required for the continuum limit) decreases the energetic cost associated with the area of the minimal 
spanning surface [23] . As the number of worldsheets of a given area grows exponentially in that area, absent 
other effects, at a sufficiently large /? the entropy of surfaces will "win out" , and the surface will delocalize 
(the position of nearby plaquettes on the surface will become uncorrelated) . This general phenomenon is 
referred to as a roughening transition |21j . 

In what follows, we show how to apply the spin foam form of dual lattice gauge theory to resolving 
each contribution to (Or„) as a local interaction between a F-connected fundamental worldsheet and a 
"background" spin foam from the vacuum ensemble. Upon doing so it is clear from inspecting the amplitude 
(and verified by numerical simulations |12j). that increasing /? also leads to rapid growth in the average spin 
on the background spin foams that the worldsheet has to interact with. Analysis of this interaction will lead 
us to conjecture that it is this growth in the average background spin and the resulting high energy cost of 
inserting large area worldsheets that provides confinement in the weak coupling limit. 

3. Existence of Worldsheet-background decompositions for SU{2) 

In this section we will start by reviewing an equivalent formulation of Yang-Mills theory on the lattice 
that uses spin foams, which are closely related to strong coupling diagrams but allow the factorization of 
the amplitude into local degrees of freedom. The reader is referred to some of the existing literature for 
derivations and details on the spin foam formulation of lattice gauge theory [i m[34l[35] and earlier forms of 
non-abelian duality P^-T in the specific case of _D = 3, G = SU{2). 

Using the spin foam formulation, we will show how a Wilson loop observable can be expressed as a sum 
over worldsheets bounded by the loop interacting with a sum over vacuum spin foam configurations. An 
explicit construction for this interaction energy can be given and will turn out to be local to the worldsheet. 

To begin, we define a spin foam for a D-dimensional lattice k and group G as follows. 
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Definition 2 (Lattice spin foam). Given a group G and a lattice k with edges E and plaquettes P, a spin 
foam is a pair of maps: a plaquette labelling P — > 7?.(G) and an edge labelling E 1{G), where 7^(G) and 
1{G) are the unitary irreducible representations of G and the intertwiners of G, respectively. 

Using the spin foam model dual to lattice Yang-Mills, the vacuum partition function can be written as 

(8) zv^Y. -^(/) = E n ^^(/(^')) n ^^(/(^)) n Avuiv)), 

f&Fv feFv peP eeE v£V 

where Fy denotes the set of lattice spin foams associated with the vacuum partition function. The amplitudes 
Ay, Ae, and Ap are functions that depend on the plaquette and edge variables that are local to each vertex, 
edge, and plaquette respectively. The set f{v) consists of the irrep variables assigned to plaquettes incident 
to V and to the intertwiner variables on edges incident to v. The set /(e) denotes the intertwiner variable 
assigned to e and the irrep variables on plaquettes incident to e by /. The set f{p) is simply the irrep value 
assigned by / to p. While the spin foam transformation may be carried out for any SU {N) , we shall restrict 
ourselves to the case of SU{2) in what follows, as methods for explicitly evaluating amplitudes and their 
asymptotics are in this case more fully worked out at present than for SU(3). 

To characterize the topological structure of those spin foams that make a non-zero contribution to Zy, we 
introduce the notions of edge and spin foam admissibility for vacuum spin foams (where no external charge 
is present). 

Definition 3 (Edge vacuum admissibility). Let an edge e G i? be given. Let n equal the number of plaquettes 
incident to e, and Vj. the vector spaces associated to each representation label ji. An edge is admissible if 
the space of intertwiners / = lnv{Vj-^ ® • ■ • ® V,„) is non-empty and the edge variable is labelled from the 
set /. A convenient graphical rule for determining edge admissibility can be given as follows. 

For every Vj- , mark dim(l/, . ) — 1 points on the boundary of a disc, with each point associated with a given 
irrep grouped beside each other. Then the dimension of the space / = Inv(V,j (g) • • • ® Vj^) is equal to the 
number of crossingless matchings in which no matching occurs between points in the same grouping, and 
each such crossingless matching^ corresponds to a different intertwiner }25| . An illustration of a crossingless 
matchings with and without U-turns is given in Figure [1] below (in this example we work in D = 3) . 




Figure 1. (a) A crossingless matching with "U-turn" (does not correspond to any inter- 
twiner) (b) Illustration of a crossingless matching corresponding to one possible intertwiner, 
given a labelling of plaquettes incident to an edge. 



Definition 4 (Spin foam vacuum admissibility). A spin foam is admissible if Ve £ E, edge admissibility is 
satisfied. 

We now turn to the modification of these conditions in the presence of a Wilson loop. We seek a spin 
foam expansion of the numerator of equation (|4]) , which we shall define to be the externally charged partition 
function Zr- First we note that away from T, the spin foam amplitude is unchanged. One consequence of 



A matching is a pairing of the boundary points of a disc by lines that connect them through the disc; the matching is 
crossingless if the hnes can be drawn in such a way that none of them intersect. 
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this is the modification of vertex ampHtudes along T. The spin foam amplitude for the charged partition 
function is thus 

feFr f&FrP&P eeEr t>GVr 



where 

(9) Avi{jU{^}v 



Milh 



_\Av{.fiv)) forw^r 
[A*y{f{v),T) forwer 

UEifie)) fore^r 
[A*E{f{e)) foreer 

For edges and vertices that are part of F, modified amplitudes and Ay are necessary (see for example [SHU] 
for details) and in the case of A*{V) depend on whether F passes straight through v or makes a turn. 

As with the vacuum diagrams, we seek to characterize the non-zero contributions to Z-p in terms of local 
admissibility conditions. As the insertion of the trace around F only introduces representation matrices 
depending on for e G F, admissibility conditions on edges away from F are unaffected. Thus for defining 
admissibility with regard to Zp it suffices to consider e £ F, as follows: 

Definition 5 (Edge F- admissibility). Let an edge e belonging to a Wilson loop F be given, let n be equal 
to the number of plaquettes incident to e, and let Vj. denote the vector spaces associated to each plaquette 
irrep label ji. We say that e is T-edge admissible if and only if the space Inv(Vi ® Vj-^ 05 Vj^ ■ ■ ■ ® Vj„) is 
non-empty and the edge is labelled by an intcrtwiner belonging to this non-empty set. 

This leads to a charged edge admissibility condition corresponding to a crossingless matching structure 
with the incident representations appearing as in the vacuum case, plus a single additional point, which we 
shall call the charged point. The condition is illustrated graphically in Figure [2] below. 

We note in passing that this definition can be applied to a loop F with multiple connected components; 
e.g. a Wilson loop observable with two disconnected loops separated by some distance, such as is used in 
computing the two-point expectation value (See Appendix [B|. 

Since the group integrals at edges away from F are unchanged, the admissibility conditions there are 
equivalent to those in the vacuum. Conjuncting the appropriate admissibility condition at each edge of the 
lattice we have the following definition for spin foam F- admissibility: 

Definition 6 (Spin foam F- admissibility). Given a Wilson loop F, a spin foam is F-admissible if and only 
Ve G E, edges belonging to F are F-admissible according to the Definition [51 and edges not belonging to F 
are vacuum admissible according to Definition [3] 

The above gives a local definition of spin foams contributing to the expectation value of a Wilson loop 
observable. The topological structure of the allowed configurations can be derived by considering (starting 
from the spin foam with all trivial labels), first all the ways in which F-admissibility can be satisfied locally 
at F, and then considering the remaining choice^. To arrive at a F-admissible spin foam, we first satisfy 
admissibility for each edge on F by specifying a matching there. Depending on how the charged point is 
matched, this will in turn determine a set of plaquettes connected to F. Denote this set of plaquettes by 
p(o), Por each plaquette in p("\ for each of its edges not attached to F there is at least one point that has 
to be routed, as the spin on the P*^") plaquettes is greater than or equal to j = ^. The situation for a given 
edge at this stage of construction is illustrated in Figure [2] below. Upon choosing a routing for each edge, 
there is another set P'^^ of plaquettes on the other end of the routing. Repeating the above construction 
will lead to a set of plaquettes that grows until all routings close, in which case P'^^^ U P'-^-' • • ■ will define a 
closed surface of (at least) half-integer charge. Note that the growing surface can fold back (not immediately 
at an edge, but after first passing through an edge) on itself; in this case the spins required to define the 
surface may exceed j = \ and the surface is said to be self-intersecting. 

We note that during any step of this construction (and after the surface is closed), an arbitrary labelling of 
plaquettes that by themselves satisfy vacuum admissibility conditions at every edge can be added (addition 



As admissibility constraints are additive, this order could be reversed without effecting the conclusion. 
"^Or more precisely, self-intersecting at a plaquette rather than a vertex or edge; this is the most useful definition for what 
will follow. 



j(p*)>0 




j(p*) 



j(p*)>0 




j(p*) > 



Wilson charge must get routed 
to some plaquette to p* 



j=l/2 



Other edges ofp* must pass non-zero 
charge to neighboring plaquettes 



Figure 2. Construction and matching conditions for P*^"^ (three dimensional example) 



of irrep labels at plaquettes) to whatever spins are present by the procedure given above. Moreover, since the 
plaquette labelling is edge admissible everywhere, there is at least one and in general multiple edge labellings 
that constitute admissible spin foams. Thus, the most general solution to the admissibility constraints for 
a Wilson loop insertion has the form of a direct sum of a spin foam in the form of a closed (possibly self- 
intersecting) worldsheet of charge bounding F, and a "background" spin foam satisfying vacuum boundary 
conditions. 

We hasten to add that while a specific choice of steps in the above construction gives a unique spin foam, 
the converse does not hold. That is, there does not exist a bijection p : Fr ^ Fy x Sy where Sr is the set of 
closed surfaces connected to F and Fy is the set of vacuum admissible spin foams defined above. In practical 
terms, this means a double sum over F-connected worldsheets and over background spin foams would result 
in an overcounting as the same spin foam in Fy would appear multiple times. Further discussion of this 
degeneracy is provided in Appendix [Al 

Notwithstanding, there are several ways to proceed. First, one can simply choose a specific decomposition 
p : Fr ^ Fv x Sr, and for each worldsheet in Sr sum over some restricted set of Fy (or vice versa). 
Unfortunately, it may be difficult to systematically organize such a sum. However, some work that is 
particularly relevant to this approach is described by Conrady and Khavkine in |16j . where each spin foam 
is associated with a unique configuration of fundamental worldsheets that are self-avoiding on a thickened 
two-complex associated with the original lattice. This approach appears natural topologically and would be 
interesting to develop further. 

Another alternative that would allow a direct sum over F-connccted worldsheets and vacuum spin foams 
is to divide each term by a degeneracy factor equal to the total number of distinct ways the associated spin 
foam can be realized by Fy x Sr- 

Finally, a (possibly crude) over-estimate for the Wilson loop expectation value could be formed by ignoring 
degeneracy when summing over worldsheets and background spin foams. Because of the possibility of 
negative amplitudes in the spin foam model, an absolute value would have to be taken to get a meaningful 
bound. Unfortunately, this over-counting of states combined with taking the absolute value (eliminating any 
interference effects) may result in too crude an upper bound to provide a useful confinement result. 

The above derivation and considerations are combined in the following statement: 

Theorem 3.1. Let a function p : Fr ~~* Fy x Sr be given. Let S^{f) be the set of s G Sr such that (/, s) is 
in the image of the map p. Let D^^{f,s) be the degeneracy factor associated to the worldsheet-background 
pair {f,s), as defined in Appendix [Al Let Sr(A) be the set of worldsheets with area equal to A in lattice 



units. Then the following relations hold: 
(10) {Or J ^leFrMf) 



E/gi.,.(E.gsg(/)J(/,^)) Mf) 

Ef^F,. \EA>A^.^ir)EsesAA)D-HLs)IiLs)] A{f) 



EfeF.Mf) 



wh- 



ere 



is the interaction factor that measures the change in amplitude associated with inserting a worldsheet s into 
a background spin foam f to obtain a new spin foam f. The numerators of ill]) are evaluated on the spin 
foam f = f + s, and can be thought of as functions of the original spin foam variables displaced by the 
presence of the worldsheet s. The functions f{p), /(e) and f{v) pick out the labels assigned by the spin foam 
to a plaquette, edge or vertex, respectively. 

The first equality is simply the spin foam expression for Wilson loop expectation values as a ratio of 
partition functions built from F-admissible and vacuum spin foams. 

In the second equality, a sum is performed over all background spin foams. For each background spin 
foam, a subset 5p(/) of all F-attached worldsheets are summed over; in general this is a proper subset 
of Sr to avoid overcounting. In the third equality, an unrestricted double sum over all background spin 
foams and all worldsheets is performed, with an inverse degeneracy factor D~^{f, s) included to correct the 
overcounting. In the final line, we have organized the sum over F-terminated worldsheets by their surface 
area (in plaquette units). The sum starts at ^Min, the area of the smallest closed surface that terminates on 
F. 

We should emphasize several interesting properties of the interaction factor /(/, s) that appears above. 
First, we observe that (for a given worldsheet s) it is an observable on the vacuum ensemble. Although 
division occurs in its definition, on the space of admissible vacuum spin foams the observable is always finite. 
Moreover, as is manifest from its definition it factors into a product of local terms, and in particular the 
plaquette term can be readily found in closed form. With the above in hand, we now turn to a possible 
confinement scenario. 

4. Confinement at Weak Coupling: A Scenario 

We describe here a general scenario for confinement using the above reformulation of the Wilson loop 
criteria in terms of worldsheets interacting with background spin foams from the vacuum ensemble. To 
begin, we introduce the notion of confinement relative to a fixed "background" spin foam from Fy^ the set of 
vacuum spin foams. A first important step is to distinguish the plaquette and vertex parts of the interaction 
factor: 

(12) ivu. s) = n AAfiv)) ' ^-(^(^)' - n Ap(/(,)) ■ 

We assume the edge amplitude can be absorbed into the vertex amplitude, which can be generically achieved 
by including a square root of the edge amplitude on either vertex of a given edge; with this convention the 
interaction term factors as /(/, s) — Iv{f, s)/p(/, s). From Theorem 13. 1[ we have 



(13) (Or J 



E/eFv Eyi>AMin(r) EsGSr(A) ^ ^{f,s)I{f,s) A{f) 



^A>AMi„(r) Z^sGSr(A) 



EfeF.m 
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where we have used the resummation by area that will be particularly well-suited for comparing effective 
tensions in what follows. 

We shall next present a sequence of arguments for how the above expression may lead to confining behavior 
in the weak coupling (/3 — > oo) limit. Our first step is to provide for the existence of (Or„} at a fixed T 
by identifying the conditions under which it is bounded from above by a convergent exponential series in 
the worldsheet area A. We then compare the leading order terms for T of different size (minimal area) and 
propose a condition under which they are bounded by an exponential decay in minimal area. 

4.1. Area damping from plaquette interaction. An elementary analysis of the plaquette interaction 
factor Ip for the SU{2) case (see [TT] for details — note we are using whole rather than half integers to label 
irreps) with the heat kernel action [29| provides a plaquette interaction that has the form 



(14) Ip{f,s) = (^e'^'"i^pesSpi'^)iP^ (^e-m-'i:pesSp{s)iSAs)+2 



^ (,-rp(J.s)A(s)^ 

where the sum is over every plaquette in the worldsheet s and we define the plaquette tension as Tp(/, s) = 
'°^^A(s('^"'"' '^here A{s) is the lattice area of s and Sp{s) is the irrep associated to the plaquette by the 
worldsheet. Note the second bracket factor in p4| goes to unity as /3 ^ oo; consequently, this factor does 
not play the crucial role that it does in the strong coupling limit. Given that the third factor can be made 
arbitrarily large on any background on surfaces of high enough self-intersection, how can a positive tension 
be provided? We see by inspection that if the jp (on the backgrounds that dominate Zv at any given (3) 
grows at the same rate or faster than /? as /3 ^ oo, than the first factor will be less than unity (provide 
tension to each worldsheet) on those backgrounds. Moreover, it will also damp the third factor (exponential 
decay in S vs linear growth in S) for sheets of arbitrarily high self-intersection. To continue our analysis to 
the rest of the interaction factor we substitute our exponentiated form of the plaquette interaction p4p into 
the expectation loop expression (I13p above: 

E,eF,. \EA>A..,.,r)E.^s.,MD-Hj.s)e-''<'^-''-'Iv(J.s)] M!) 



For a given /, an infinite number of above unity interaction terms (the summands appearing in square 
brackets) in (|15p above would ruin convergence, or at best give conditional convergence and associated 
ambiguities. Thus, we would like to have below unity interaction factors on all but a finite set of sheets (at 
least on those background / that aren't suppressed in the weak coupling limit). 

Inspection of (jlSp reveals further conditions: Tp(/, s) at large (3 must not only be positive, it must be 
sufficient to overcome the growth in the number of surfaces of a given area, which also grows exponentially 
(is bounded from above) in the area [12] with some associated negative tensior0 , i.e. as e^aiA-A^i^) ^ 
note in passing that Tq is dimension dependent and so in principle may lead to different tensions or even 
departures from confining behavior depending on dimension. 

As the degeneracy factor is always less than unity, the remaining tension after accounting for the growth 
in surface number e^°^ in each area class has to be sufficient to overcome any above unity interaction factors 
from the vertex amplitude. To isolate the tension coming from the plaquette amplitude, we consider a related 
model where /y = 1 (note that Ay in A{f) is not set to unity) and make the following definition: 

Definition 7 (Area damping for plaquette interaction). A plaquette interaction (and the plaquette ampli- 
tude from which it is derived) is said to be area damping if 



(16) 



< oo. 



As both vertex and plaquette interactions factor over the area of a worldsheet, what is gained by analyzing 
them separately? In the G = SU{2), D — 3 case it was found in [IT] that Iv{f,s) 1 for large spin 
backgrounds (characteristic of the weak coupling limit) and self-avoiding worldsheets. Consequently Ip 



^Although we suppress it in the notation To = To{A) carries some A dependence. 
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rather than ly would have to provide confinement for this class of worldsheets and backgrounds. The large 
(3 form of Ip (jl4[) provides such a means to provide effective tension, assuming the typical spin grows fast 
enough with increasing /3 as mentioned above. This motivates our conjecture that for Yang-Mills theories in 
general gauge groups and dimensions, it is the Ip that provides an effective tension as described above. If 
this conjecture is true, than area damping is a necessary condition for the existence of the expectation value. 
Given the area damping condition, the behavior of the vertex interaction becomes decisive; we turn to this 
next. 

4.2. Stability of vertex interaction. For the area damping from plaquette interaction to provide confine- 
ment within our proposed scenario, the behavior of the vertex amplitude that would lead to negative overall 
tension from (above unity interaction factors) must be bounded from above to a degree determined by the 
other tensions (on all but a finite set of sheets). We start by making the following definition: 

(17) /y(/,s)^e(/,s)7'(/,s)e-^"(/-^)-4(^), 

where e(/, s) is the overall sign of vertex interaction factor, the quantity T-'{f, s) is the perimeter factoi^ 
defined as V{f,s) = | Yiver -^^if^ ^)\ ^^"^ defined implicitly by the above. Unlike the plaquette factor, 
we note /v'(/, s) can in principle introduce a negative sign, as it involves spin network evaluations that are 
not strictly positive. 

We next introduce effective area-averaged tension f(/, s) as 

^ ^ ^ ~ (,ToiA-AMin) f.To{A-AMi,-,) 

To compare the this effective tension from the interaction factor with the "negative" tension arising from 
growth in surface number, we consider the following form of the expectation value: 

(19) {Or J = ^^^"^ ^ ^ ^ 



which suggests the inequality 

(20) f{f,A)-To>0 

as a necessary condition for a contribution to the sum over A to be below unity (for a given background 
spin foam /). If this inequality is satisfied for all A for a foam / (and the infinum of r over all the A is 
non-zero), then the vacuum amplitude A{f) is multiplied by a sum of terms having the form of a constant 
of order e~^0'4Min^ g^^,^ ^ factor exponentially decaying at a rate proportional to the area. This leads to the 
following definition: 

Definition 8 (Vertex Stability). The set of ordered pairs of {f,A) where inequality (|'20p holds is said to 
exhibit vertex stability. 

It may be the case that the inequality is violated on some finite subset of areas — in which case the sum of 
interaction factors still converges. Additionally, if there are any / £ Fy for which there are areas for which 
inequality violation occurs, they should make vanishing contribution to the expectation value in the large 
/3 limit. If this were not the case it may be very difficult (although not impossible) to organize the leading 
order of the sums to find an area law. An example of suppression of such ill-behaved world-sheet background 
sectors can in fact be shown for the "sub-asymptotic" configurations (essentially those containing low spin) 
in the D = 3, G = SU{2) model [H]- 

4.3. Balance conditions and a scenario for confinement. The weak coupling limit (Or„) presented 
in (fT9|) subject to the inequality ([20|) is highly suggestive of a convergent interaction factor and a possible 
area law. That is, if A(f) were strictly non-negative one could form an upper bound with the form of a 
series of decaying exponentials in the area, and compare the leading order contribution associated to T of 
increasing size. Unfortunately, positivity does not hold for A{f) in a general spin foam model for Yang-Mills. 



We factor out a perimeter factor as the vertex amplitudes have a different form along T, as discussed above; this is often 
of practical utility when working with a specific model. 
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Hence, before proceeding further we must deal with the non-positivity of the ampHtude. One approach to 
account for this is to write the ansatz (fT9l) above as follows: 



(21) 



(Or J 



-To Ami n 



A>AMi„ ' 



,(To^tU\A))A 



(E/eF+ - E 



-To AmIii 



A>AMin ' 



,{To-f{f.A))A 



\A{f)\ 



\feF- 



Ef^F.Mf) 



E 



+ e 



-ToAf. 



A>AMin ' 



,{To-fif,A))A 



\A{f)\ 



E/,F.|^(/)I 



'E 



-To A.Min 



E 



E/eF- l^(/)l 



(To-t(/,A))A|^(-j^|\ 



■ 



where we have used the fact that E/Gi="v ^('^) E/gf^ ^r(/) are absolutely convergenlQ and can thus 
be broken into two sums E/s-F^ "^^^^ = E/eF+ ^ E/eF- 

We next make the following conjecture: the interaction factor for any given worldsheet s, averaged over 
negative amplitude states Fy approaches its average over positive amplitude states Fy as (3 oo. That is, 



(22) 



'E 



-ToAt. 



E 



„(To-r(/,A))A 



^(/)r 



'E/ 



-To Ami n 



E/.F,.|^(/)I 

E 



A>A. 



,(To-t(/,A))A 



as P ^ oo. 



As an immediate consequence we then have from 
(23) {Or J ^ + 

which allows us to write 



Ef,F^\Aif)-Ef,F,m)\ \ 

EfeF^lMn-EfeF-lMDl) 



(24) 



(Or 



E 



E/eF,tl^(/)l 



'To Am in 



E 



< 



A>AMin 



-t'A 



E 



l^(/)l 



A>AMia 



where r* = infj^T»+ ^ [''■(/, A) — Tq] is the infinum (greatest lower bound) of the effective tension over areas 

and Fy , the set of positive amplitude vacuum spin that give positive tension; it is assumed to exist and 
is non-zerc|l. The term E(P) represents contributions which may have zero or negative tension but go to 
zero as /3 ^ cxj. An example of such are sub-asymptotic contributions discussed in [11]. In a theory where 
confinement is true, we refer to Fy ^ as the dominant sector. 

A final step in showing confinement is to compare the leading order terms in (Or,„ ) for different loops T 
of increasing size (lattice area of minimal surface). A useful construct is to consider the series of background 



This can be shown for any given (3 due to upper bounds on the vacuum and F vertex amplitudes. 

Otherwise terms of arbitrarily small tension would be present in the continuum limit; the theory would not be confining. 
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averaged expectation values: 



(25) «,(A) = ±^^^^v 



EfeF^ \A{f)\ 
in terms of which we can write 

(26) (Or„) = e-^'"^™" ^ ar{A). 

We next need a way of comparing ar(^Min) for T of increasing size. Suppose we make a generalization of the 
infinum r* used in ([24|) . As r* in that case is defined for a specific F, we write T*(r) and take the infinum 
over some set of loops T of interest, G- 

(27) feEEinfe[T*(r)]. 
If f > 0, then the following relation will be true: 

(28) an (^Mi„(r)) < e-^s^™" for aU TeG. 

Note that the condition pS)) by itself permits all the manner of fluctuations in (Or„) as F is increased, as 
long as they all take place within a decaying envelope with fall-off given by f . It is however not difficult to 
show monotonic decrease in the Q;r(^), both for a fixed T and varying A and for varying T. Such a result 
puts additional bounds on the extent of fluctuations in the f. 



4.4. Critical Remarks. The area damping and vertex stability conditions represent non-trivial statements 
regarding the spin foam form of the gauge theory in the weak coupling limit. For any given model, further 
analysis will be required to examine the validity of these propositions. Notwithstanding, if these proposals 
prove false the alternatives for realizing confinement are quite constrained. We briefly mention a few of the 
possibilities here. 

It may be the case that for certain gauge groups, the area damping condition fails, i.e. the plaquette 
amplitude provides insufficient tension relative to the surface growth factor on backgrounds that dominate 
the /? — > oo limit. If this "vertex damping" scenario is true, it would be quite remarkable as it would place 
very specific constraints on displacements of vertex amplitudes with arbitrarily large spin arguments, whose 
contribution to Zy grows as f3 oo. 

A second scenario more closely related to the one proposed is that the plaquette and vertex amplitudes 
are by themselves insufficient but combined can overcome the surface number growth. Because neither type 
of interaction factor is alone sufficient, one could term this a "cooperative" confinement scenario. 

In considering the final steps of the above, the balance condition calls for particular scrutiny, as it was 
motivated largely for its convenient mathematical properties. If it proves to be false, a weakened form may 
still be true. In this weakened form, the balance condition may hold true only for a subset of spin foams, in 
which case the remaining spin foams would have to be shown to be suppressed or exhibit area law behavior 
by some other line of argument. If such balance conditions prove false (or impractical to show), it may 
be fruitful to consider other means of passing from the form (|19p to an area law inequality, such as by 
employing a stationary phase analysis. Such an analysis may be particularly appropriate in the asymptotic 
setting where analytic formulas for vertex amplitudes are available, as in [TT] . 

The presence of the degeneracy factor in defining the effective tension is also somewhat unsatisfying, as 
we at present don't have a good estimate for it. Thus, we don't a priori know that configurations with large 
vertex interactions which might otherwise spoil confinement aren't damped by the degeneracy factor. A 
particularly well-chosen decomposition p (see Theorem 3.1 for definition) might offer a better understanding 
and is the topic of ongoing work by the author. Notwithstanding, if sufficiently strong bounds on the behavior 
of vertex interactions can be found, they may be sufficient (combined with estimates on other interaction 
factors) to provide at least an area law bound. Such a bound likely would not estimate the actual tension 
very well, due to a lack of degeneracy factor estimate (other than that it is less than unity). 

A final comment relates to our treatment of the perimeter factors. Our analysis does not rule out a 
non-trivial perimeter law, as meeting the conditions of area damping, vertex stability, and the existence of 
the infinum over F yields an area-law bound (rather than equality). 
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5. Conclusions 



Based on some elementary considerations that allow Wilson loop expectation values to be expressed as 
a double sum over worldsheets and vacuum spin foam backgrounds (modulo a degeneracy factor), we have 
attempted to map out some of the properties we expect to hold at weak coupling in order to provide an 
area-law decay. Specifically, the concepts introduced in this paper can summarized as follows: 

(1) Wilson loop expectation values can be expressed as a certain type of observable having the form of 
a sum over worldsheets bounded by F. Each worldsheet contributes an interaction factor that is a 
function of variables local to the worldsheet. 

(2) The interaction factors into a plaquette interaction and a vertex interaction. It is conjectured that 
the large /? form of the plaquette interaction provides the critical amount of tension required for area 
law behavior. 

(3) An inequality relating the effective tension provided by interaction factor and the worldsheet number 
in each area class is conjectured to hold on vacuum spin foams that dominate the large (3 limit. The 
contribution to the expectation value of those spin foams and area classes where the inequality is 
violated must be shown to go to zero as [3 ^ oo. 

Put in more descriptive terms, we propose that a necessary feature for confinement at weak coupling is 
that the typical vacuum configurations at large (3 exhibit sufficiently high spin relative to (3 such that the 
interaction factor (a product of both plaquette and vertex interaction factors) averaged over worldsheets of 
a given area provides an effective tension greater than the "negative tension" provided by the "entropy" 
(number of) lattice worldsheets bounding F of that area. 

Qualitatively, this mechanism is suggestive of the dual superconductor model in which the Yang-Mills 
vacuum "expels" electric flux. It may be possible to develop this connection further — spin foams can be 
thought of as path integrals over spin networks, which provide a complete basis for the physical Hilbert space 
of the lattice regularized gauge theory [7l[8] . Intriguingly, in the spin network basis, the "electric" part of 
the lattice Hamiltonian is diagonal [T'j'F and the eigenvalues are products of Casimirs in the irrep labels, 
which appear explicitly in the heat kernel form of the plaquette amplitude which leads to the plaquette 
interaction factor. We expect that, roughly speaking, the plaquette amplitude arises from the electric part 
of the lattice action, however further work on the relation between Hamiltonian and Lagrangian formalisms 
is likely required to make this statement more precise. 

At this point we would like to emphasize some differences with the confinement at strong coupling. In going 
from strong to weak coupling the confinement mechanism shifts from a high tension worldsheet embedded in 
sparse (low-spin) backgrounds, to that of a low tension worldsheet that is forced (again, by virtue of having 
to end on F) to interact with vacuum spin foam where high spins form a percolating background of spin 
which grows with increasing (3. In the weak coupling limit, partially sub-asymptotic configurations (that can 
exhibit negative tension through the vertex interactions) become exponentially suppressed, and thus don't 
spoil the confinement. 

Due to the limit nature of these arguments, which take place in the opposite extremes of strong and 
weak coupling, an intermediate region of "roughening" is not ruled out in regions of the phase diagram (e.g. 
coupling and lattice size) in which neither strong or weak mechanisms are sufficient to provide a substantial 
tension. Through the use of dual spin foam algorithms, it should be possible to quantify the relative influence 
of these two mechanisms at a given coupling. 

To further assess the validity of this proposed confinement mechanism, a number of additional results will 
be needed, the details of which depend on the space-time dimension and gauge group considered. While 
the case of SU{2) in three space-time dimensions taken up in [TT] seems so far consistent with the proposed 
picture, the case of SU{2) on a four-dimensional lattice remains to be considered in more detail. A formula 
for the four-dimensional vertex amplitude has recently been given in [13|; being built up of a sum of 6j 
and 3j symbols it may be possible to find a suitable asymptotic description simply by substitution, or else 
through a more direct approach. In both three and four dimensions, more careful control of the asymptotic 
expression will be required to handle the case of self-intersecting worldsheets interacting with asymptotic 
spin foams. 

To extend the analysis to SU{3), suitable asymptotic formulas for this case will be required in order to 
show that the vertex interaction trivializes or at least provides some bounded amount of negative tension. To 
the author's knowledge, asymptotic formulae analogous to Ponzano-Regge asymptotics for SU{2) 6j symbols 
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Figure 3. Multiple decompositions of the same F- admissible spin foam (top) into world- 
sheet (middle) plus background (bottom). 



have not yet been found. In general, the computation of 6j and 3j symbols and their recoupling theory is 
considerably less developed than the SU (2) case at this time. However, we hope that some of the work to 
date [6 u24n25j will provide a strong starting point for eventually understanding the asymptotics of SU{'i) 
spin networks. 

While the above problems may present considerable technical challenges, we believe that the main result 
of our work can be stated clearly: assuming that Yang-Mills theory truly describes our world and that 
confinement is true, the weak coupling behavior of vertex interactions and plaquette interactions relative to 
the growth rate in surface number with area are subject to non-trivial constraints. These constraints suggest 
very specific lines of inquiry, some of which are identified in here and in llj. 

In pursuing a deeper understanding of these constraints, we are encouraged by some powerful new meth- 
ods and results involving the asymptotics of 6j symbols and spin networks, from both the mathematics 
literature [Tl[Tni[17] and in the context of the spin foam program for quantum gravity jl8ji20„26j, for which 
the asymptotics of spin foam amplitudes also play a crucial role. 

Acknowledgement. The author would like to thank Dan Christensen, Florian Conrady, and Igor Khavkine 
for valuable discussions that have influenced this work. 

A. Degeneracy of Spin Foam Decompositions 

As discussed above, it is particularly natural to view contributions to (Or„ ) as composed of a worldsheet 
bounding F together with some "background" spin foam associated with the undisturbed Yang-Mills vacuum 
(no external charge). A summation over all worldsheets and all backgrounds will (quite seriously in most 
cases) over-count the underlying configurations Fr- 

A simple example of this degeneracy is the spin foam shown in Figure [S] At the top is illustrate a F- 
admissible spin foam. It can be decomposed as a worldsheet connected to F (dotted line) which intersects a 
closed surface belonging to Fy ■ Alternatively, it can be decomposed as the direct sum of a worldsheet raised 
upwards in a way that conforms with the boundary of the box and a smaller closed surface. Finally, the 
configuration where the worldsheet is pulled downwards and a smaller closed surface is fit over top is also 
equivalent. 

In the above case, a F-admissible spin foam had three separate decompositions. This leads us to the 
following definition: 
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Figure 4. Example of a pyramid-like F-admissible spin foam leading to a recursive set of degeneracies. 



Definition 9 (Degeneracy factor). Let 7 be a F-admissible spin foam. As discussed in Section 3, there 
exists at least one ordered pair of vacuum spin foams / and s such that 7 = f + s. Define the degeneracy 
factor to be the number of distinct such decompositions for a given 7. 

We note in passing the £'(7) can also be written as D{f, s) = D{f + s) and because it is non-zero the 
function D^^{f, s) is always well-defined. 

In theory, one may attempt to compute a degeneracy factor by finding how many closed background spin 
foams are sliced "in half" by the worldsheet. The pyramid-like construct shown in Figure H] illustrates how 
this may be difficult, as additional degeneracies can be generated from an initial slice. 

B. Mass Gap and Generalized Wilson loop criteria 
We give here the following (standard) definition of the mass gap observable: 
Definition 10 (Mass gap observable). 

(29) (Onr^w) ^ (OriOr^) - (Or^) (Or.) 

where Fi and F2 are the boundaries of two parallel plaquettes (one above the other) separated by distance 
X. The one-point expectation values (Opi) and (Or.) defined as above with F^ = F^. The product 
expectation value (OrjOr.) is defined as follows: 

/n n \ - ^rir. _ / {UeeEdge) Tr (Heen ^M) Tr (Heer, DM) e-^.^-^^^'') 
(dU) (CnOr.) = 



Zv J(n^^^d5,)e-2:,.P5(,,) 

As in the case of confinement (where the Wilson loop consisted of a single connected component), an 
analysis of admissibility at the two components Fi and F2 allows contributions ^rir. to be resolved into a 
background spin foam upon which are superimposed one or more surfaces bounding the F. Specifically, we 
can write a sum of the form given in Theorem 13. II 



(31) (Or.r.(.)) = ^^^^^ Aif) 

^ i:,e..(E...g^,^,^,(,)/(/,^))^(/) 

where the only difference with the case of a single connected Wilson loop component is that the set S'r has 
to be suitably generalized to account for the disconnected topology of F. One finds that S'r for a two-point 
function consists of the following types: 

(1) An (arbitrary) worldsheet with both Fi and F2 as boundary components. 

(2) Two worldsheets that bound Fi and F2. 



As in the case of confinement, the worldsheets may be self-intersecting. In the second case, the two separate 
worldshccts may also overlap. Note that worldsheets of Type 2 will become insensitive to separation in the 
continuum limit; the non-trivial behavior of interest thus lies in the Type 1 worldsheets. 

We note in passing that by all the same arguments provided in Section 3 and 4, terms contributing 
to {Orir2{x)} ■will have the canonical form of background amplitudes multipled by a sum of worldsheet 
interaction factors that provide an effective tension. The main difference is that surface "entropy", or 
number of Type 1 worldsheets of a given area may have a different growth rate (Tq in the notation used 
above) on account of their different topology With this in mind, an entirely analogous argument to that of 
Section 4 can be made with the conclusion that if the inequality ([^0]) is satisfied (with the possible exception 
of a vanishing weight subset in the continuum limit), one obtains a converging series in Type 1 worldsheets of 
increasing area, the leading order term of which will be a Type 1 worldsheet corresponding straight cylinder 
connecting the two plaquettes. Eventually we hope that similar analyses can be applied to arbitrary n-point 
functions, possibly leading to a rigorous definition of the quantum field theory. 
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